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Abstract 


For a finite mean supercriticial Bellman-Harris process, let Zt be the num¬ 
ber of particles at time t. There exist numbers (the Seneta constants) 
such that XtZt converges almost surely to a non-degenerate limit. Fur¬ 
thermore, xt oc e“ l3t C(e ~ f3t ), where f3 is the Malthusian parameter, and 
C is slowly varying at zero. We obtain a characterisation of the slowly 
varying part of the Seneta constants under the assumption that the life¬ 
time distribution of particles is strongly non-lattice. 
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1 Introduction 


We consider a supercritical Bellman-Harris process {Z t }t> o with offspring dis¬ 
tribution {7Tfc}£L 0 and life-time distribution G. In words, G(t ) is the probability 
that a newborn particle survives at least until time t, and nk is the probabil¬ 
ity that once it splits into a number Z + of progeny, it will split into exactly 
k of these. We denote by f(s) the corresponding generating function (PGF) 
/(s) := E(e sZ +) = YlkLo n kS k , and set 


h(s) 


1 ~f(s) 
1 — s 


( 1 ) 


Let y, := /i(l) be finite. It is known that there exist ‘constants’ \t (the Seneta 
constants) such that, on the set of non-extinction, \tZ t converges almost surely 
to a non-degenerate random variable Z. With Ft{s) the PGF of the distribution 
of particle numbers at time t. and F-t(s) its inverse, an immediate candidate 
for the Xt’s is 

Xt ■■=■ Xt($) := - log F_ t (<r) (2) 

for some s e (q, 1), where q is such that h(q) = 1 (such a q exists and is 
unique because of supercriticality). This is because the Laplace transform of 
the random variable Z (with y as the dummy variable) is 

R(y) :=: R,(y) := lim F t (e~^ y ) , (3) 

t —KX> 


which by definition of Xt equals c for y = 1. Hence Z is non-trivial in the sense 
that its Laplace transform is neither 0 nor 1. To get a feeling for how quickly 
the Xt’s tend to zero, recall that Z t grows essentially as e^ 4 as t —> oo, where f3 
is the Malthusian parameter, that is, 



1 

I* 


It is therefore natural to conjecture that Xt = e~ l3t C(e~ ,3t ) for some slowly 
varying (in e~ ,4t ) function C. Under this assumption, it is easy to derive an 
equation for the Laplace transform R of the random variable Z\ Since F t fulfills 
the integral equation JT] 


F t (s)= (1-G(t))a+ f foF t _ u (s)dG(u), 

Jo 

it follows immediately that 

R(y) = lim (l — G(t))e~ XtV + [ f o F t - u (e~ Xt ~ u ^ Xt ^ Xt ~ u ^ v ) dG(t) 
t~> °° J o 


foR(ye-^)dG(t) 


( 4 ) 
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because of dominated convergence and the fact that Xt/Xt-u = e~ /3u except 
for a factor which tends to 1 as t —> oo. Furthermore, since X* e_/3t is certainly 
very close to zero for t large enough, we may approximate 1 — R(e~ l3t ) ~ 1 — 
F t (e~ Xte Pt ) by E(Z t )xte~ ,3t ~ Xt , so that 1 — i?(e _/3t ) would also be a natural 
guess at the value of the Seneta constants. (By ~ we mean that the ratio of 
the two quantities is bounded from above and away from zero.) This guess is 
indeed a good one, as has been established by Schuh [|J for the case of G being 
non-lattice, and makes it natural for us to consider 

/>00 

X(t) := e 0t (1 - R{e~ pt )) = ho R(e~^ t+u ^)X(t + u)e~ 0u dG(u) . (5) 

Jo 

Our idea is to assume h, R , and G as given and treat Equation as an equation 
in the unknown function X. It will turn out that under the assumption that G 
is strongly non-lattice, which is to say that 


lim inf 

1 0 1 —►oo 


poo 

1- / dG(t) 

Jo 


> 0 , 


we will need no more than elementary renewal theory to derive the following 

Theorem. Suppose the life-time distribution G of particles is strongly non¬ 
lattice, and that E(Z+ logZ+) = oo. Then, with 

f‘00 

v ~ 1 te~ l3t dG(t) , (6) 

Jo 

we have 

log [ (pt - h{l - e~ 0u )) du , (7) 

h Jo 

where ~ means that the ratio of both sides tends to 1 as t —+ oo. 

The gist of the theorem is that both sides of 0 are asymptotically equivalent 
with a constant of proportionality equal to 1 instead of only 0(1). Note that X 
does not depend anymore on g on this scale. We also observe that under what 
might be called Uchiyama’s (K condition: 


p - h(l - s) = (-logs) “£(— logs) , 


where a > 0 and C is slowly varying at infinity, the theorem reduces to 
log X(t)~- f {I3u)~ a £{/3u) du ~ -- V t 1 ~ a C(pt), 

h Jo (! - u)P a P 

which, if we set G(t) = S T (t) (the Dirac mass at some splitting time r), (3 = 
log p/r, and v = r _1 , is what Uchiyama’s Theorem [B] claims for the ordinary 
Galton-Watson process. (Time should then be measured in units of r.) Let us 
now turn to the 
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2 Proof of the Theorem 

To begin, write 

1 f 00 

X(t) = - ho R(e- 0u )X(u) dGp{u - t) , 
h Ju=t 

with 

Gp(t) :=ii[ e-P u dG(u), 

Jo 

and introduce the ‘renewal function’ 

oo 

U 0 (t) :=£<#(*), 


( 8 ) 

(9) 

( 10 ) 


where GV is for the Tfold convolution of Gp with itself. (The standard definition 
would be to also include a Dirac mass at zero.) By renewal theory [5], we have 
for t > 0 , 

1 + Up(t) = ut + Up{t) , (11) 

where Up(t) >0, v is as given in the theorem, and 


dUp(t) = — t 2 dGp(t) . 


( 12 ) 


^ poo pvAt 

h 
l 


M J v—0 J u= 0 

pt pV 


Since G /3 has an exponential tail, and is strongly non-lattice together with G, 
Up has an exponential (right) tail as well p|- Consider now 

[ X(u) dUp(u) = ~ [ [ hoR{ e -P v )X(v)dGp{v-u)dUp(u) 

J u =0 M J u—0 J v=u 

h o R(e~ 0v )X{v) dGp(v - u) dUp(u) 

h o R(e~ 0v )X(v) dGp{v - u) dU 0 (u ) 

h o R(e~ 0v )X(v) dGp(v - u) dUp(u) 

h o R(e-P v )X(v) d(Gp * Up){v) 

1 f°° 

- ho R(e~ 0v )X(v) d(Gp * Up(v) - Gp * Up(v - t)) 

M J V=t 

h o R(e~ 0v )X{v) dGp(v - u) d(Up(u) - Up(u - t)) , 


M J v—0 J u—0 

dr r 

M Jv=tJu—0 

1 


poo pV 


M J v—0 
+ ■ 

1 

M Jv=tJu=t 
since, by Equation mr 

Up(u + t) = vt + Up(u ) + Up{u + t) - Up(u) . 
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But 


Gp * Up(t) = Up * Gp(t) = Up(t) - Gp(t) 
by definition of the renewal function, so we obtain 

pt i pt 


' u—0 


1 I 

X{u)dUp(u) = — / h o R{e~P v )X(y) dUp(v) 

M J v =0 

hoR(e~ 0v )X(v) dGp(v) 
h o R(e-P v )X(v) d(Up(v) - Up(v - t)) 

P OO 1 pOO 

/ hoR(e~ pv )X{v)dGp(v) + - / h o R( e - 0v )X(v) dGp(v - t) 

J V=t M «/ V=£ 

1 /»00 /‘•y 

- / / h o R(e~ Pv )X{v) dGp(v - u) d(Up{u) - Up(u - t)) . 

d J v=t j u=t 

The second and forth term on the right-hand side of this equation add to — X (0), 
by Equation 111 111 . The fifth term, by the same equation, is simply X(t). As for 
the remaining term, we have 


M J v=0 

+ i r 

d Jv=t 
2 /*°° 

M Jv= 

p OO PV 


1 


p OO pV 


M JV=t JU=t 

p OO -I POO 


h o R(e~ 0v )X(v ) - it) d(Up{u) - Up(u - t)) 

POO -i POO 

/ - hoR{e~ pv )X{v)dGp{v-u)d(Up{u)-Up{u-t)) 

Ju=t M Jv=u 

X(u) d(Up(u) - Up(u - t)) , 


so that 

f 

Ju =0 


X(u) dUp(u) = 


+ X(t) - X(0) - 


1 r 

- / h o R(e~P u )X(u) dUp(u) 
M Ju= 0 

/-0 ° n-ho R(e- 0U ) 


J U=t M 

If now we make use of o once more, we finally obtain 

r t 


X(u) d(Up(u) — Up(u — t)) . 


where 


-[ (p-hoR(e~ l3u ))X(u)du = X(t)-X(0) + X(t)-X(0), (13) 

M Ju =0 

3 p — h o R[e~P u ) 


With 


r{t) := 


X(t) := 

m 

X(t) 


■ X(u) dUp{u — t) . 


(14) 


u—0 


H-ho R{e~^ t+U '>) X(t + it) 
p X{t) 


dUp(u ) > 0 , (15) 


6 







we can write down a ‘solution’ for this as 




(i-o 


1 + g(0) 




l-i — h o R(e 0U ) 
1 + a(u) 



( 16 ) 


We can assume that under E(Z+ log Z+) = oo, the integral in the exponent 
diverges (see, for instance, Corollary 3.1 in [3], or Lemma 2 in [Sj). Therefore, 


log X(t) - ~ I 
d Jo 


H — h o R(e 0U ) 
1 + cr(u ) 


du 


certainly if cr is bounded for all sufficiently large t. We show that this is the 
case: Since X{t) is slowly varying in e~ 0t , it follows by Potter’s Theorem 2] 
that for arbitrary A > 1, <5 > 0, there exists T\ := S) such that 


X(t + u) 


< Ae S0u 


for all t>T\. Thus we obtain for all such t , from Equation GSt 


a(t) < A 


fj, — h o R(e 0t ) 
d 



dUp{u) , 


(17) 


which is of order R(e~ 0t ) for suitably chosen <5, because Up has an 

exponential tail (<5 < 1 already suffices). In particular, a(t) —> 0 as t —> oo, 
which readily implies that even 


\ogX{t) ~ — / (p, — h o R(e~ 0u )) du =: y(t) . (18) 

d Jo 

Now, because of 03 , and due to the divergence of the integral in m, there 
exists r 2 such that X(t) > 1 for all t > r 2 , and 


y{t)-y(T 2 )>- f (n~h o(l-e 0u ))du, 
M J T2 


by monotonicity of h. But y(j 2 ) is finite, hence 

y(t) 


lim inf-j- 

t - > °° vp- 1 f Q (fi — h o (1 — e ~ 0u )) du 


> 1. 


(19) 


On the other hand, because X(t) is slowly varying in e 0t , we can choose A > 1, 
S > 0, and 73 < 00 such that Ae~ S0T3 X(t 3 } < 1, and 

y{t) - y(r 3 ) < J (m - h o (! - e -/3(1-,5) “)) du 

for all t > 73 , by monotonicity of h , and Potter’s Theorem again. But 3 ^( 73 ) is 
finite, hence 


lim sup 


y(t) 


< 


1 


t^oo vp, 1 f* (// — h o (1 — e /3u )) du 1 — 5 


( 20 ) 
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which together with © and the definition of y concludes the proof of the 
theorem, since <5 is arbitrary. □ 
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